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ABSTRACT 

New generation ground and space-based CMB experiments have ushered in discoveries of 
massive galaxy clusters via the Sunyaev-Zeldovich (SZ) effect, providing a new window for 
studying cluster astrophysics and cosmology. Many of the newly discovered, SZ-selected clus- 
ters contain hot intracluster plasma (IcT e > lOkeV) and exhibit disturbed morphology, indica- 
tive of frequent mergers with large peculiar velocity (v > 1000 kms -1 ). It is well-known that 
for the interpretation of the SZ signal from hot, moving galaxy clusters, relativistic correc- 
tions must be taken into account, and in this work, we present a fast and accurate method for 
computing these effects. Our approach is based on an alternative derivation of the Boltzmann 
collision term which provides new physical insight into the sources of different kinematic 
corrections in the scattering problem. This allows us to obtain a clean separation of kinematic 
and scattering terms which differs from previous works. We also briefly mention additional 
complications connected with kinematic effects that should be considered when interpreting 
future SZ data for individual clusters. One of the main outcomes of this work is SZpack, a 
numerical library which allows very fast and precise (< 0.001% at frequencies hv < 20kT y ) 
computation of the SZ signals up to high electron temperature (kT e ^ 25 keV) and large pe- 
culiar velocity (v/c ^ 0.01). The accuracy is well beyond the current and future precision of 
SZ observations and practically eliminates uncertainties related to more expensive numerical 
evaluation of the Boltzmann collision term. Our new approach should therefore be useful for 
analyzing future high-resolution, multi-frequency SZ observations as well as computing the 
predicted SZ effect signals from numerical simulations. 
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1 INTRODUCTION 

Free electrons residing inside the deep potential wells of galaxy 
clusters scatter photons of the cosmic microwave background 
(CMB), causing a spectral distortion which is commonly referred 
to as the Sunyaev-Zeldovich (SZ) effect. The thermal/random mo- 
tions of electrons in the hot cluster atmosphere s lead to the thermal 
SZ (thSZ) effect l lZeldovich & Sunvaevll 19691) . exhibiting a y-type 
spectral shape related to the up-scattering of CMB photons. The 
bulk/directed motion of the electrons, on the other hand, causes a 
temperature shift in th e direction of the cluster, kno wn as the kine- 
matic SZ (kSZ) effect l lSunvaev & Zeldovich|[T980l) . 

The SZ effect has long been realized as a powerful tool 
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to l earn about the f ormation of struc t ures in the Universe 
(see iRephaelil Il995al ; iBirkinshawl 1 19991 ; ICarlstrom et"ail |2002|, 
for reviews). Since the effect is independent of redshift, the 
SZ signal is ideally suited to discover massive galaxy clus- 
ters out to high redshifts. New generations of ground-based 
(e.g. iMuchovei et alj200l;IVanderlinde et al.l20ld:lMarriage et alj 
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I201I and space-based JPlanck Collaboration et al 
experiments have lead to discoveries of several 
selected clusters. Current cluster samples have already pro 
vided unique constraints on the dark energy equati on of 



vided uni que constraints onjtne dark energy equation ol 
state (e.g.. Wang et alj |2004|; iHutsil l200fj| ; Ivikhlinin et alj 12009] : 
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20111). the Hubble c o nstant (e.g., 
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l2002l ; lBonamente et alKOOfj). and non-Ga ussianity in the primor- 
dial matter density field dFolevetal .11201 ll) . 
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Many of the newly discovered SZ-selecte d clusters contain 
very h ot intracluster plasma (kT e > lOkeV) dWilliamson et alj 
1201 ll : iKorngut et all l201lh and exhibit disturbed morphol- 
ogy, indicative of frequ ent mergers with large peculiar veloc- 
ity (v > 1000 kms~') faalverson et"ai] 120091 : iMa et alj |2012| ; 
iMroczkowski et al]|2012l) . The 'Bullet' cluster is one well-studied 
examples of a high veloc it y merging sy s tem, w i th v ^ 4000 kms -1 
1 Markevitch et alj |2002|; iMarkevitch] 120061: ISpringel & Farrar l 



20071: Milosavlievic et alj 120071 : iMastropietro & Burkerj 1200^ : 
Lee & KomatsitolCh . 'El Gordo' is another spectacular example 
of a massive, merging cluster at z = 0.87 fille d with very high 
tempe rature gas up to kT e - 22keV ± 6keV jMenanteau et alj 
|2012|) . Furthermore, IMroczkowski et all d2012l) recently measured 
the SZ effect for a triple merger system, finding shock heated 
gas with temperature AT e > 20keV and large peculiar motion 
(v ^ 3000 kms~'). Because of the high electron temperature and 
large kinematic terms, interpretation of the observed SZ signal for 
such systems clearly must take into account relativistic corrections. 

These corrections have bee n investigated theoretically by 
several independent groups (e.g.. Iwrightj|l979l: lRephaelilll995bl: 



Challinor & Lasenbvl Il998l lltoh et al.l Il998l: ISazonov & Sunvaeyl 
19981 : iNozawa et all 1 19981: iMolnar & Birkinshawl I1999T 



Enfflin & Raised I2OO0I: iDolgov et alj l200ll: [Shimon & Rephaelil 



20041 : IChluba et alj 120051 : IPoutanen & Vurml hoiol) . and various 
analytic approximations already provide a good description of the 
SZ signals in different ranges of temperatures and frequencies. 
However, a number of high-resolution SZ experiments, including 
ALMAQ, CARM^E CCAli, and MUSTAN(fJ are underway or 
planned, promising a dramatic increase in sensitivities, spatial 
resolution, and spectral coverage over the next few years. These 
upcoming SZ experiments should enable a host of new mea- 
surements of important cluster properties, in cluding the electron 
temperature of the intraclu s ter medium (ICM) ( Pointecouteau et al] 
1 19981 : lHansenetal]|2002l : IZemcov et alj|2012h , the peculiar ve- 
locity and interna l bulk and turbulent gas motions of clusters 



dNagai et all 120031: ISunvaev et all 120031 : biego et alj 120031) . and 
non-equi librium electrons produ c ed by merger and a c cretio n 
shocks dMarkevitch & Vikhlininl 120071 : iRudd & Nagail 120091) . 
to name a few. In order to properly interpret these upcoming 
measurements, it is critical to develop a fast and precise method for 
computing the SZ effect signal, ideally with accuracy well beyond 
the precision of future SZ experiments. 

In this work, we present a fast and accurate method for 
computing the relativistic corrections to the SZ effect. The ap- 
proach is based on an alternative derivation of the Compton col- 
lision terms, obtained using explicit Lorentz-transformations of the 
SZ signals from the cluster frame into the observer frame. This 
method allows us to clearly identify the source of different kine- 
matic correction terms. The cluster's peculiar velocity, /? c = v/c, 
makes the electron distribution function effectively anisotropic in 
the CMB rest frame. One way of dealing with the problem is there- 
fore to use the Lorentz-boosted Bolt z mann equation in the CMB 
rest frame (e.g.. INozawa et all 1 19981: IChallinor & Lasenbvl [l999l : 
INozawa et al .120061) . Alternatively, one can treat the problem in the 
cluster frame. Upon Lorentz-transformation of the isotropic CMB 
into this frame, the CMB spectrum becomes anisotropic (e.g., see 
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IChlubal201 with motion-induced monopole (oc /? 2 ), dipole (oc f} c ) 
and quadrupole (oc /? 2 ), due to aberration and Doppler boosting. 
These anisotropies are then scattered by hot electrons, imprinting a 
distortion on the CMB. Using the Lorentz-invariance of the photon 
occupation number and line-of-sight number of scatterings one can 
readily obtain the SZ intensity signal in the observer's frame. The 
advantage is that in the cluster frame complications related to kine- 
matic corrections (aberration, retardation and time-dilation effects) 
and the cluster geometry (see Sect. 14.5) can be avoided. 

The crucial generalization is that we consider not only 
the scattering of the radiation monopole, but also dipole and 
quadrupole scattering by hot electrons. The required kinetic equa- 
tion for the scattering of anisotropic radiation in lowest order 
of the electron temperature was recently derived by IChluba et alj 
d2012h : however, here we include higher order temperature correc- 
tion to this problem. We obtain a reformulation of the frequency- 
dependent basis functions (Sect. 13.1.11 and 14. 1 .4b that allow very 
precise (< 0.001% at frequencies hv < 20kT y ) computation of the 
SZ signals up to high electron temperature (kT s ^ 25keV) and 
large peculiar velocity (v/c ^ 0.01). Our new approach pushes 
the precision for the SZ predictions far beyond the sensitivity 
of future SZ observation, eliminating the need for more expen- 
sive numerical integrations of the Boltzmann collision term. It 
should therefore be useful when analyzing future high resolution, 
multi-frequency SZ data for individual clusters or whe n computing 
the p r edicted SZ signals fr o m numerical simu l ations ( Nagai et al] 
20031: ISunvaev et alj|2003l: iDiego et alj|2003l: iDolag et alj|2005l: 



Battaglia et alj|20ld) . 



With our new derivation, we also show that the kinematic 
corrections to the SZ intensity sign al given here differ from 
previously obtained expressions (e.g., ISazonov &Sunvaev|[i998l: 
INozawa et all 1 19981 , 120061 : IShimon & Rephaelil I2004T ). This dif- 
ference is related to the interpretation of the scattering optical 
depth integral, which was not explicitly addressed in the ear- 
lier calculations (see Sect. 14.41 . Although, the differences are 
small, they could still be relevant for instance when using the 
SZ effect to confirm the redshift-scaling of the CMB tempera- 
ture jBattistelli et alj|2002l:lHorellou et al|2005l : lLuzzi et alj2009t 
Ide Martino et alj 2012: lAvgoustidis et al.ll20i20. an ide a that was 
suggested long ago dFabbri et alj|l978t lRephaelilll980l) . We also 
briefly mention additional complications connected with kinematic 
effects on the measurement that should be taken into account when 
interpreting future SZ data for individual clusters (see Sect. |4.5l l. 
However, a detailed analysis is beyond the scope of this paper. 

One of the main outcomes of this work is SZpackB a numerical 
library which allows explicit computation of the SZ signals using 
the full Compton collision integral, but also fast approximation of 
the numerical result by means of the improved set of frequency- 
dependent basis functions obtained here. This new approach al- 
lows calculation of the SZ signals with precision < 0.001% over 
a wide range of parameters at practically no computational cost, 
overcoming limitations (see Sect. 13.11 for m ore details) of previ- 
ously derived analytic approximations (e.g . IChallinor & Lasenbvl 
1 19981 : lltoh etaill 1998! : ISazonov & Sunvaevlll9981) to describe the 
SZ signals from hot, moving clusters. With respect to other fast rep- 
resentations of the Bolt zmann collision by means of extended ta- 
bles or fitting functions 1 Nozawa et al .120001 : lltoh & Nozawal2004l: 
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IShimon & Rephaelifcoojh our approach provides an interesting al- 
ternative, with the derived basis functions being directly informed 
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by the underlying physics of the problem. We furthermore include 
the kinematic corrections in this representation, describing the SZ 
intensity signal in the cluster frame and then obtaining the final 
result by Lorentz-transformation. This ap proach also allows us to 
include the effect of the observer's motion dChluba et al.l2005l) in a 
simple manner, without extra approximations. 



2 ANISOTROPIC COMPTON SCATTERING 

In this section we present a brief derivation of the required kinetic 
equation for Compton scattering of photons by hot thermal elec- 
trons. For now, we neglect the effect of peculiar motion and only 
consider relativistic correction in different orders of the electron 
temperature, T c . However, we describe the scattering of anisotropic 
photon distributions in the rest frame of the SZ cluster. This gen- 
eralization is needed to include kinematic corrections up to second 
order in the peculiar moti on using explicit Lore ntz transformations, 
as recently pointed out bv lChluba et alj d2012l) . 



2.1 The Boltzmann collision term for a resting cluster 

The time evolution of the photon phase space density n{t, v, y) 
at frequency v in the directior0 y under Compton scattering 
can be described using the Boltzman n equation (compare with 
iBuchler & Yuehlll97d : lftoh et al.ll 19981) 



dn(t, v, y) 
dt 



(i) 



where d 2 y' is the solid angle spanned by the incoming photon, T 
is the statistical factor, and dcr/ AO! denotes the differential scat- 
tering cross sectiorQ. We assumed that the incoming photon field is 
spatially homogeneous, which for the problem of interest here is a 
valid approximation as long as the number of scatterings is small. 

To compute the SZ effect one can neglect the effect of electron 
recoil, since the energy of CMB photons, hv, is much smaller than 
the energy of the fast electrons inside the cluster, i.e., hv <k ym c c 2 . 
Here y = 1 / - [S 2 is the Lorentz factor and ft the dimensionless 
velocity of the moving electron. In that limit the Compt on scatter- 
ing cross section reads (e.g., see Ijauch & Rohrlichll 19761) 



dcr 

dfi' 



3<x x 



1 

y-K 



v y-K 1 1 \ v y 2 K 2 



(2) 



with k = 1 -ftp, a sc = 1 - p sc , where p sc = y' ■ y is the cosine 
of the scattering angle between the incoming and outgoing photon. 
Furthermore, p = • f and p! = ■ y are the direction cosines 
of the angle between the scattering electron and the incoming and 
outgoing photon, respectively. The Thomson scattering cross sec- 
tion is denoted by cr x as 6.65 x 10~ 25 cm 2 . The frequency ratio of 
the incoming and outgoing photon is given by 



l-/3p 



(3) 



Since the temperature of the electron gas obeys kT c <K m e c 2 , 
Fermi-blocking is negligible. In addition stimulated scattering can 
be omitted so that the statistical factor T may be approximated as 



In the following bold font denotes 3-dimensional vectors and an addi- 
tional hat means that it is normalized to unity. 

7 More concisely, the factor related to the M0ller relative velocity was ab- 
sorbed in the definition of dcr/ dfl' . 



f « f n' — fn, where the abbreviations / = /(p), /' = /(p') 
for the electron, and n = n{y,y), n' = n(v',y') for the photon 
phase space densities were introduced. The electron phase space 
density is isotropic and may be described by a relativistic Maxwell- 
Boltzmann distribution, 

^ (P) " 47r(m c c) 3 £ 2 (l/#e)#e ' (4) 
where A^(l/£! e ) is the modified Bessel function of second kind with 
dimensionless electron temperature 6 C = kT^jm^c 2 , N c is the elec- 
tron number density, and e(p) = -\J 1 + ?/(p) 2 = y denotes the di- 
mensionless energy of the electrons with the dimensionless mo- 
mentum 77(p) = |p|/m c c = y/3. 

With the above definitions, the statistical factor can be cast 
into the form T ~ f[e Axc n' - n], where we defined x c = M- and 
Aj c = x' e — x c . Introducing x = hv/kT y and A v = (v' - v)/v, the 
statistical factor can be rewritten as 



rim - r) + Xir ** 5 * n(x ' ^ 



(5) 



where we set Ax c = (r r /T c ) x A v x 0, neglecting terms that are 
multiplied by the CMB to electron temperature ratio, T r /T s , which 
for typical clusters is of order ^ 10~ 8 . 

We now define the moments of the energy shifts 

1 C dcr „. . M 



r 



r dcr 

J dfi 7 



m-^YUy')d 2 y'd 3 p 
k\ 



(6) 



/V* c cr T 

over the scattering cross section. Here Yi,„(y') denote spherical har- 
monic functions. Since Yi(-,„)(y') = (-1)"' Yf m (y'), it directly fol- 
lows /jf_ m > = (-l) m (^,)*, so that only the integrals with m > 
ever have to be explicitly computed. For small # c these moments 
are all frequency-independent polynomials, which can be explic- 
itly computed up to some order, k max , in 8 e . As we also see below, 
If = Ai Yi„,(y) by symmetry, where ?k is a temperature-dependent 
function. In Appendix [C] we give the results for the moments of 
interest up to 0(S)} ) for monopole scattering, and (9(6!*) for dipole 
and quadrupole scattering 

If we write n(x,y') as a spherical harmonic expansion with 
spherical harmonic coefficients n lm (x), then the required Boltzmann 
equation takes the simple form 

^»-n(,,« + ZZ4^, (7) 

*=0 l,m 

where we introduced dr = cN e cr T dr. In the optically thin limit, this 
expression is formally valid for any anisotropic, yet spatially uni- 
form radiation field that is scattered by hot electrons and in which 
electron recoil and stimulated scatterings can be neglected. For the 
temperature correction to the SZ effect one only needs to consider 
the scattering of the monopole part of the radiation field, but when 
considering kinematic corrections, anisotropies in the photon dis- 
tribution must also be included. 

With modern computers it is straightforward to carry out the 
full 5-dimensional collision integral, Eq. Q 3, num e rically . Here we 
use the quadrature rules according to IPattersonl (1968) for each 
of the independent 1 -dimensional integrals. These rules are fully 
nested and converge rapidly with only a few refinements. Below 
we simplify the collision term by reducing the number of integrals 
and by leveraging symmetries of the problem in different situa- 
tions. To confirm the accuracy of the results we always compare 
to those obtained with the full 5-dimensional collision term. With 
the full collision term it is also possible to compute the SZ signals 
for more general incoming radiation fields with anisotropies. The 
corresponding routines can be found in SZpack. 
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3 THERMAL SZ EFFECT FOR HOT ELECTRONS 

To compute the temperature correction to the SZ effect for t <k 1 
we have to insert n ss ;i P i = l/(e A - 1) into the Boltzmann equa- 
tion, Eq. Q. Up to n th order in f9 e we need to compute the deriva- 
tives x*3*Wpi for k < In + 2. These derivatives can be expressed in 
several ways, some of which are summarized in Appendix [A] For 
explicit computations of the relativistic corrections it is enough to 
directly use the expressions for x*3*rc P i as they are, taking care of 
the order of different terms to improve the numerical stability (see 
Appendix [A). Together with the moments in Table IbTI they deter- 
mine all correction terms up to the desired order in temperature, in 
our case (9(f?"). The spectral distortion caused by the thSZ effect in 
the single scattering approximation is then given by 



0.2 
0.15 

0.1 
0.05 


-0.05 
-0.1 
-0.15, 



e e = o.oi 
e c = 0.02 
e e = 0.03 

Y approximation (10th order) 




10 20 30 



A«, h (x,y)« J]y (A) K,(x) 



(8a) 



cr T cdt « f3f +1 Ar = f^> <0) , (8b) 



where we introduced the generalized y-parameter, y (k \ and the 
line-of-sight optical deptfQ At = J N e cr^cdt. For the second 
and third equalities of Eq. l !8bt , we also assumed that the clus- 
ter is isothermal along the line-of-sight. Alternatively, one can 
consider the distortion caused by a single scattering event in one 
small volume element of the cluster and then sum each contribu- 
tion along the line-of-sight. The functions Y„ are determined by 
Y* = Z£f ajt°**^Hoo C*) y 00 with «oo(x)Too = "piOO, and the coef- 
ficients a? summarized in Table IB ll 

With the definitions given in Appendix {A), in first order of 
the electron temperature we have 



An 



At 



e e [4xd x + x 2 ^ 2 .] « P1 = e e g Lccoth 1 1 j - 4] = ej , 



with Q(x) = xe x l[e x — l] 2 = —xd x n n ( x). Th is is the well-known 
result given bv IZeldovich & Sunvaevl (fl96^ for the thSZ effect, 
which can be directly derived usin g the Kompaneets equation 
dKompaneets!ll956l : IWevmand[l965l) . In second order of the elec- 
tron temperature, we find 



An 



12) 



At 



10xd t + — x 1 d\ + -x 3 d\ + —x 4 dt 
2 5 1 10 1 



tpi = %Y X 



This expression is equivalent to Eq. (2.27) inlltoh et al.l dl998h and 
the last term in Eq. (12) of ISazonov & Sunvaevl dl998l)~ as can 
be easily shown using the expressions, Eq. dAlt . We confirmed 
higher order ter ms by comparing with the expressions given by 
lltoh et alj l fl998l) . Our expressions for the temperature corrections 
to the thSZ effect alone are als o equivalent to those formulae given 
by|sh imon & Rephaelil J2004I) . although here we obtained two ad- 
ditional orders in the electron temperature. We also derived tem- 
perature corrections up to 20 lh order in electron temperature, but 
improvements in the convergence were only marginal. 



Unless noted otherwise henceforth we will refer to 'optical depth' as 
quantity that is defined inside the cluster frame. With this convention the 
optical depth is directly characterizing the number of scatterings along the 
photon's world line, and hence cleanly separates the effect of scattering and 
kinematic terms on the SZ signal, as we explain in detail in Sect. 14.2.11 



Comparison with numerical result 




6 -0.01 (10th order) x 100 




8 e = 0.02 (10th order) 




... 8 -0.03 (10th order) 
j 

/ i 




' [i 

Crossover frequency j 
i 


' ; i 1 

i ' ; l ' 
■1 Ml 

■ i i r 
i • ; 1 L 
; l it 1 

■ i \ ^. 



Figure 1. Thermal SZ effect for an isothermal electron gas with different 
temperatures, f? c = (kT^lm^c 1 ). In the upper panel we present the change 
in the intensity, in dimensionless form, A// At k jc 3 An/Ar = x 3 Y(x), where 
Y(x) = Ztfn Yk(x), neglecting terms 0(8]}). The lower panel shows the 
relative difference with respect to the full numerical result. For e = 0.01 
we multiplied it by 100, demonstrating that the SZ signal is represented to 
better than =: 0.01% precision over the shown frequency range. However, 
for clusters with kT e > lOkeV the Taylor series expansion with the basis 
functions Ft does not perform well at x > 10. Notice that the spike around 
the crossover frequency, x CI ~ 3.83, is because the SZ signal vanishes there. 
The numerical integrals were carried out using SZpack. 



3.1 Improved representation of the spectral distortions 
caused by the thermal SZ effect of hot clusters 

Although the moments /L have previously been computed to very 
high orders in temperature, it is kn own that the associated ap- 
proximation co n verges very slowly dChallinor & Lasenbvl Il99l: 
lltoh et alj 1 19981 : ISazonov & Sunvaevl Il998l) . This is because for 
high electr on temperature (> lOkeV) the width of the scattering 
kernel (see lSazonov & Sunvaevll2000l for properties of the Comp- 
ton kernel) becomes sufficiently large that the absolute frequency 
shift (Ax ^ x V2(? e - x/5) caused by Doppler broadening exceeds 
unity at x > 1/V2#T - 5. Although the relative change in fre- 
quency is still small, the simple Taylor series expansion of npi(x') 
around x starts to converge slowly: even when taking up to 10 th or- 
der correction in f5 c , the formula provide a rather poor description 
of the SZ signal at high frequencie^] However, the frequency range 
x - 10 - 20 is very relevant for interpreting observational data. 



9 We found that pushing to much higher order in f> e (we computed the 20 th 
order) the match degrades significantly, (cf. Fig.[T). 
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It is therefore important to compute the thSZ effect from hot clus- 
ters precisely even at these high frequencies. Furthermore, higher 
order temperature corrections are significant near the crossover fre- 
quency, x cr =! 3.83, and improved formulae would be desired. 

One straightforward possibility is direct numerical integration 
of the Boltzmann collision term, Eq. QJ. Since for the thSZ effect of 
hot clust ers the number of integrals can be analytically reduced to 
two(e.g..lWrightll979l : lEnfilin & Kaiseill2000l : iDolgov et alj200ll ; 



iNozawa et al. 20091) . the line-of-sight numerical integration for a 
single temperature is very fast (a few seconds on a standard laptop) 
and a full analytic description in terms of a Taylor series can be 
avoided. However, when analyzing upcoming high resolution SZ 
data along many lines-of-sight, the inversion of the SZ signal to 
temperature and electron profile requires a large number of inte- 
grals to be taken. In that case it is beneficial to have a simpler and 
faster method for computing the SZ signal. 

One possibility is to produce a simple 2-dimen sional table in 
(6 e , x ) and then interpolate or fit the results (see e.g., Nozawa et alj 
120001 : lltoh & Nozaw3l2004k IShimon & Rephaelill2004h . but for an 
accurate representation of the collision term (< 0.1% say) one 
needs a rather fine grid, rendering this approach suboptimal . In par- 
ticular, it might be useful to look for a different approach when 
kinematic corrections are also included. 

Here we follow an alternative approach in the middle of the 
two extremes: one can use a Taylor series in 9 C for the Boltzmann 
collision term, but leave rcpi(x') as it is, carrying out the remaining 
integrals numerically. This procedure defines a set of frequency- 
dependent basis functions, Z k (x), that do not suffer from the same 
limitation as the Y k (x) expansion, while still being informed by the 
underlying physics of the scattering problem. 



3.1.1 Computing the improved basis for the thermal SZ effect 

To compute this new set of basis functions, we first reduce the di- 
mensionality of the collision integral by realizing that for isotropic 
incoming radiation the azimuthal integrals can be carried out an- 
alytically (see Appendix |B1 for details). The electron temperature 
enters the problem via the electron distribution function /(p), de- 
fined in Eq. (4]l. If we factorize /(p) as 



/(P) 



{ K 2 (l/e c )8j 47r(m e c) 3 ' 
with Ay = 7 - 1 , we can define the basis function!"'! 



dV b*e 



dp dp' d k 6 e 



■ [npi(x') - npi(x)] dp dp'ifdr] (9) 



4/r 



£ 2 (l/6y£> c (27r6y 3 ' 2 



i-H £ 
2 



345 
— ( 

128 



+ 0(6l) 



which can be computed numerically for some fixed temperature, 
# e = 6 c . The azimuthally averaged cross section d 2 cr / dp dp' is 
given by Eq. ( IB lab . Then the distortion caused by the thSZ effect 



10 One could in principle carry out either the integral over p or p' in addi- 
tion by defini ng the logarithmic frequency shift s = ln[(l-/?/j)/(l-/?//)] (cf. 



IWrightll 19791) . but the resulting expressions are more complicated and nu- 
merical integration over 3 dimensions remains sufficiently straightforward 
with high precision so that we stop here. 




Figure 2. Comparison of the basis functions Y k and Z k for 6 e fl = 0.03. We 
scaled both sets by appropriate factors of k and ft c to make them comparable 
in amplitude. One can clearly see that the new set of functions, Z k , exhibits 
fewer oscillations as k increases. This leads to more stable convergence for 
higher orders in ft, . 



of a hot cluster is determined by 
A/Wx^^Az^Oe), 



(10) 



A;' 



1/,) 



(6> e - 6» e ,o) dT * 



We) 



(0 e -0e,o) AT, 



Me e .o) Af(0 e ,o) 

where in the last step we assumed that the electron gas is isother- 
mal over the distance dl = cdt. This assumption is always justified 
when choosing dl sufficiently small, i.e., by considering the scat- 
tering effect for a small volume element of the cluster. The deriva- 
tives d k e e^ hy ^ c can be performed analytically for any k. In Fig. [2] 
we compare the first few functions Z k with Y k . One can see that the 
new set of functions, Z k , exhibits fewer oscillations as k increases. 
This leads to more stable convergence for higher orders in S . 

In Fig. [3] we illustrate the performance of the new expan- 
sion. We computed the basis functions, Z k , at reference tempera- 
ture (? c _o = 0.03 up to the 20 th derivative in # e using SZpack. Clearly 
the full numerical result is represented very accurately (with preci- 
sion comparable to 0.001%) over a wide range of frequencies and 
temperatures. We found that the expansion works well even at tem- 
peratures of 9 e > 0.05. At low temperatures, C < 0.01, it is better to 
rely on the asymptotic expansion based on Y k , since errors start to 
become larger. However, the combined description with Y k and the 
new functions Z k allows very precise computation of the thSZ effect 
in an extremely efficient way. A combination of both approxima- 
tions is implemented in SZpack. Alternatively, one could produce 
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Figure 3. Same as Fig. \T\ but for slightly higher electron temperature and 
using the expansion, Eq. j 10j . for different orders. Clearly the convergence 
is excellent at all frequencies of interest. 



a set of basis functions Zj around different reference temperatures 
and in this way compress the representation of the full integral in an 
economic way, while achieving even higher precision at the same 
time. This is especially important when considering effects in sec- 
ond order of the cluster's peculiar motion (velocity of = 3000km/s 
say) for which relative precision better than ^ 10~ 4 is required. 



4 KINEMATIC CORRECTIONS TO THE SZ SIGNAL 

To derive the modifications caused by the motion of the cluster 
with respect to the CMB and the observer, we describe the scat- 
tering process in the rest frame of the cluster, where the incom- 
ing radiation field is anisotropic due to kinematic effects. We can 
use the formalism developed in Sect.|2]to compute the introduced 
distortion by thinking of the problem as scattering of monopolar 
through quadrupolar part of the photon distribution. Transforma- 
tion of the SZ signal into the observer frame can then be achieved 
by simple Lorentz transformation (Sect. |4~2l and 14. 3t , without ad- 
ditional loss of precision. We restrict the main derivation to kine- 
matic corrections up to second order in the clusters velocity, j8 c , 
but neglect the effect caused by terms in 0(01) and scattering of 
primordial large-scale CMB anisotropics. We briefly discuss these 
effects in Sect. 14. 1.31 arguing that close to the crossover frequency 
these only contribute at the level of 0.1% relative to the kSZ sig- 
nal. Inside the cluster frame, the computation of the SZ signal is 
straightforward, and complications related to aberration, retarda- 
tion and time-dilation effects, as well as geometrical aspects can be 
avoided. These, however, enter the problem when transforming to 
a general observer, as we discuss in Sect. 14.51 



4.1 SZ signal in the moving cluster frame 

As a first step, we transform the photon distribution from the CMB 
rest frame into the cluster frame which is moving with a peculiar 
velocity j8 c . Aligning the z-axis with the dire ction of the cluster 's 
motion, it is straightforward to show that (cf. JChluba et al.|[2012h 



Pi 

"Zoo (*c) * noo + — 



"00, 



«l (Xc) ~ Q-l/Jc^c^WoO 



(11) 



in the cluster frame. In the chosen coordinates the photon distri- 
bution is azimuthally symmetric about C , with all other spherical 
harmonic coefficients vanishing in 0(01). We neglected any intrin- 
sic large-scale anisotropies of the CMB and set «oo = V4tt n P \(x c ). 
Also, (X\ = 1/ v3 and a 2 = 1/(3 VB"), and x c = hv c /kT y , where v c is 
the photon frequency evaluated in the cluster frame and T y defines 
the CMB monopole temperature in the CMB rest fram 

As Eq. i ll It shows, in first order of fi c a dipolar anisotropy 
is induced by the cluster's motion relative to the CMB rest frame. 
The spectrum for this part of the photon distribution corresponds 
to a temperature shift, ff(x c ) = -x c <9 Vc n P1 = x c e* / (e* - l) 2 . In 
second order of /3 C a quadrupolar anisotropy and a small correc- 
tion to the monopole appears, w here both Doppler boosting and 
the aberration effects contribute ( Challinor & van Leeuwenl 120021; 
Kosows kv & Kahnias hvili 1 20 111 ; lAmendola et all 1201 ll ; IChlubal 
201ll ; IChluba et alj|2012» . The spectrum of the quadrupolar part 



is Q(Xc) = x\Si ftpi = FoC^c) + 4^(jc c ), which has a y-type depen- 
dence. Similarly, the correction to the monopole exhibits a y-type 
distortion. Higher order correction terms (i.e., the motion-induced 
octupole, oc 0\, and hexa-decupole, oc p*) can be easily obtained, 
but for typical peculiar velocities (0 C ^ 10~ 3 - 10~ 2 ) these are very 
small (see Sect. l4,l,5l l. Also, the associated Taylor series in /? c con- 
verges quickly, which is in stark contrast to the expansion in small 
electron temperature and frequency shift discussed above, which 
only converges asymptotically. 

Inserting Eq. ( lilt into Eq. Q and collecting terms we have 



drf(x c ,f) 
<9t c 



"I 



inn -V^v 



1 + yh 



3 c * 



"pi 



+ 0£ 



AP w xA c n PI + I k w x^d^xAjip, 



(12) 



AP 20 x 2 c dln PI + £ 4 xl^xldlnn 



Here P 2 (x) = [3X 2 - l]/2 is the second Legendre polynomial and 
Hi is the direction cosine of the angle between j8 G and the outgoing 
photon, f c , evaluated in the cluster frame. Furthermore, we defined 
A/JJ, = I? - 1, and the moments = it/ Yi m (f c ), which are only 
functions of the electron temperature, f^' = kT^/m c c 2 , fixed in the 
cluster frame. Also, dr = N^crjc dt c defines the Thomson optical 
depth, with electron number density, A^:, again in the cluster frame. 

Equation dl2t is still the general expression, describing the 
scattering of CMB photons in the cluster frame, written in a com- 
pact way. The required moments are summarized in Appendix [C] 
Any observer that moves with a velocity j8 c relative to the CMB 
measures a spectral distortion in the direction of the cluster that 
is given by this expression. By including as many moments, /* , 



The value of T y can be determined, a point that we clarify below. 
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as is required to converge for a given cluster temperature, this is 
formally the exact result, provided that the value for the peculiar 
motion of the cluster allows truncation of the series in /J c at second 
order. If higher orders in /3 C are needed then additional moments TL 
with / > 2 must be computed. However, the series in /3 C converges 
rapidly, so that Eq. d!2t should be enough. As expected, when trun- 
cating the moment-hierarchy at some finite temperature, we again 
encounter problems of convergence for hot clusters. We will show 
how to overcome these limitations below. 

We also mention that, as expected, the SZ signal is indepen- 
dent of the azimuthal angle, <f> c , of /? c with respect to the line-of- 
sight, reflecting the symmetry of the scattering process. For the 
distortion only the projections of the anisotropic photon distribu- 
tion onto Yio matters, once the z-axis is re-aligned with the direction 
of y c . This statement is no longer true when the polarization signa- 
ture caused by the kSZ effect is considered JSunvaev & Zeldovichl 
Il98d :ls azonov & Sunvaev|[l99^ : in this case <* c defines the direc- 
tion perpendicular to the linear polarization plane. Internal motions 
of the ICM can cause interesting effects in the polarization signal, 
in pri nciple allowing to study the transverse intra-cl uster velocity 
field dChluba & Mannheim! [20021 ; iDiego et alj|2003l) . However, a 
treatment of the polarization effects, including temperature correc- 
tions, is beyond the scope of this paper. 



4.1.1 Effect on the monopole spectrum inside the cluster frame 

The first sum of Eq. d!2t describes the effect of scattering on the 
monopole seen in the cluster frame, with second order kinematic 
corrections to the Planckian spectrum included. The first contribu- 
tion to this sum is the thSZ effect, which affects the purely Planck- 
ian part as discussed in the previous section. The other two terms 
oc /3 2 can be included as corrections to the functions, Y k . Since we 
have a simple and numerically stable method for computing deriva- 
tives of the Planck spectrum (see Appendix IA2t . we can obtain 
the associated spectral function, 

marized in Table IB ll With the simple operator relations 

x k d k x xd x = kx t d x + x* + 1 d k ; 1 (13) 

x^d 2 = k(k - l)x*aj + 2fcc* fl 3* ri + x* +2 5f 2 , 

the effect of scattering on the monopole part of the photon distribu- 
tion, n c = iiC )0 (x c )Yoo, is given by 



Yf' n (x) using the moments sum- 



-^2(^)*[n(x c )+/j 2 y, ki "(x c )] 



(14) 



: ^ 

— ^ ,-. , 


A 




,1 I I i 
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. 2n+2 
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Figure 4. First few functions l*™, Z) km and gj™. We rescaled by appropri- 
ate factors of k and 8 e = 0.03 to make them comparable in amplitude. 



The coefficients are given in Table IB ll This expression defines 
the correction of the SZ signal in second order of y8 c arising from 
the scattering of the monopole part of the photon distribution in the 
cluster rest frame. With the definitions of Appendix lAll we can, for 
example, write 

kin 50 14 11 1 - 

Yf* =^-&+-rY -—xSm + *) + -xSX{Ug + x c ) (15) 
J 3 o o 

for the lowest order term in electron temperature. Here we intro- 
duced X = xcoth(f). The first few functions Y^ a are shown in 
Fig. [4] As for the Y k -basis of the thSZ effect we can observe strong 
oscillations for larger values of k, indicating that the associated se- 
ries in 9 C converges slowly. The collection from the motion-induced 
monopole is the smallest among all kinematic corrections to the SZ 
signal, as we discuss below. 



4.1.2 Effect on the dipole and quadrupole part of the spectrum 
inside the cluster frame 

The next two terms in Eq. l !12t are caused by scattering of the dipo- 
lar and quadrupolar anisotropy in the cluster frame. The contribu- 
tions oc A/ ; q include the effect of Thomson scattering (no energy 
exchange between electrons and photons) and temperature cor- 
rections to the Compton scattering cross section. Comparing with 
Eq. I ll It it is clear that these terms leave the spectrum of the in- 
coming radiation unchanged, but only modify their amplitude by 
scattering photons in and out of the line-of-sight. According to 
the moments given in the Tables lB2l and |B3l in second order of 



the electron temperature we have A/ 1 



Al° 

^20 



1 



mr and 



9/10 - 1 8% + 4§(0e) 2 - An observer in the rest frame of 
the moving cluster sees a dipole and quadrupole spectrum with 
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n\ (x c ) and rf 20 (x c ) given by Eq. i ll It . When looking towards the 
cluster the photon fluxes in the dipolar and quadrupolar parts of 
the photon distribution are reduced in both cases, since scatter- 
ing out of the line-of-sight dominates. For the dipolar anisotropy 
in the cluster frame the leading order term desc ribes the kSZ, 
An = ^p c yig(x c ) dSunvaev & Zeldovichl [l980h . while for the 
quadrupole scattering we have the quadratic kSZ effect (qkSZ), 
An = - ] |AT t: j e c 2 /' 2 ^)x c ^(x c )coth(x c /2) = -^At^P^)^). 
Here we introduced <3(x) = 4ff(x) + Y (x) = §(x)X(x). 

The other terms in Eq. (I12> describe changes in the spectrum 
of the dipole and quadrupole caused by the scattering process, with 
energy transfer from the electrons to the photons. For the total effect 
on the dipole part of the photon distribution we therefore have 



Kinetic SZ effect from dipole part for ]l B = 0.01 



gn"(x c .,y c ) 
3t c 



£(x c ) + ^£(^)*Z) k '"(x c ) 



(16) 



where we used the operator relations, Eq. ( 113) , and d 1 -^ according 
to Table lB2l Similarly, for the effect on the quadrupole we find 



dnii(x c ,y c ) 
<9t c 



10 



. 2n+2 



2 C 2 ]«pi, (17) 



with defined according to Table lB3l The first few functions D k: 
and Q\ m are shown in Fig. [4] We also have 



38 12 2 
■yg-y)'„ + -i c ^ + i c ) 



(18) 



So'Vc) = -^9 + -^Y - ^x c g(6g + x c ) + ^x c £l(12£ + x c ) 

in lowest order term of the electron temperature. 

Among all kinematic corrections to the SZ signal, the kSZ 
effect in first order of the cluster velocity is the largest. Temperature 
corrections to the kSZ can reach ^ 10% - 20% at high frequencies 
and are therefore larger than the correction from the qkSZ effect. 
Temperature corrections to the qkSZ effect are already very small, 
but still larger or comparable to the signal related to the motion- 
induced monopole (Sect. [47TTT}- However, as we show below, it is 
easy to include all these effects using an approach similar to the 
thSZ effect case, at no additional computational cost. 



4.1.3 Comparing with numerical results 

It is straightforward to numerically compute the spectral distor- 
tion in the rest frame of the cluster. By inserting the transformed 
occupation number, n c (x c ,y c ) = npi(y c x c (l + f) c ■ y c )), into the 
Boltzmann equation, Eq. Q, one can directly integrate the colli- 
sion term to all orders in /?<. and 6 C . We implemented this case as 
one of the routines in SZpack to confirm the results obtained with 
the simplified approach discussed below. This also allows us to in- 
clude the scattering of primordial CMB anisotropics in the cluster 
rest frame, but both higher order terms in /? c and the scattering of 
CMB anisotropics cause very small additional signals, which are 
only noticeable close to the crossover frequency, where the thSZ 
effect vanishes (see Sect. 14. 1 .5t . 

We start by considering the effect of scattering on the dipolar 




J I 
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Figure 5. Kinematic correction to the SZ signal caused by the motion- 
induced dipole (oc /3 C ) in the cluster frame for different temperatures. In 
the limit 8 C = we have the usual kSZ effect. For 8 e = 0.05 we also show 
the result of the asymptotic expansion, Eq. | |161 . with 4 th order temperature 
corrections included, indicating slow convergence at high frequencies. 



anisotropy in the cluster frame. We can ease the numerical inte- 
gration by using the symmetries of the scattering process: with re- 
spect to the outgoing photon direction y L we only have to consider 
the azimuthally symmetric part of the photon distribution function, 
ni(x',f ) = y8 c yuj;x'fl v /npi(x')yu s c, with fi x = f ■ y c , since the other 
contributions all average out. For the dipolar part of the photon dis- 
tribution we therefore have 



c^(x c , y c ) 

3t c 



=Aa£ mxc) + c km (x c ,^)\ 



C ki "(x c , (£) = -Mj e- & ^ §{x' c ) d^ d^Vdn, (19) 

where d 2 rr, / d/j dyu' is given by Eq. dB lbt . The dependence on the 
electron temperature arises from the integral C km (x c , QF), which ac- 
counts for additional redistribution of photons (up-scattering) over 
frequency by energy exchange with the thermal electrons. 

The results of the numerical integration are shown in Fig. [5] 
for different electron temperatures. The case 6 S = corresponds to 
the normal kSZ effect. Temperature corrections are visible only at 
high frequencies around x ^ 1 and 10, reaching ^ 10% - 20% rela- 
tive to the kSZ signal. For f? c = 0.05 we also show the result of the 
asymptotic expansion, Eq. l |16t , with 4 lh order temperature correc- 
tions included. At higher frequencies convergence gets slower and 
including more orders did not further improve the agreement with 
the full numerical result. On the other hand, for low temperatures 
(6> c = 0.01) we found good agreement of the expansion, Eq. d!6t . 
with the full numerical result. We also confirmed numerically that 
the full 5-dimensional collision integral for the dipole part is cor- 
rectly represented by Eq. dl9) . 

The corrections in second order of /? c can be obtained in a 
similar manner. For the motion-induced distortion of the monopole, 
defining A4(x) = Q(x) + Y (x), we find 



/3iZ km (x c ,0;) 

d2<T .-Ay/eg 



(20) 



Z km (x c , 6F) = ^ J e- & ^ [MCX) - M*c)] ^ V n 2 dn 
using Eq. iBlat and ©. Notice that Z ki "(x c , 6£) = for -> 0. For 
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Figure 6. Kinematic correction to the SZ signal in the cluster frame to sec- 
ond order of fS c for different temperatures. The upper panel shows the effect 
related to the motion-induced quadrupole. In the limit e = we obtain the 
qkSZ effect. For e = 0.05 the result of the asymptotic expansion, Eq. J17K 
with 4 th order temperature corrections included is also shown. — The lower 
panel shows the effect related to the motion-induced monopole. It vanishes 
in the limit 8 e = 0. For 8 e = 0.03 we also give the result of the asymp- 
totic expansion, Eq. 417t . with 8 th order temperature corrections included, 
indicating slow convergence at high frequencies. 



the quadrupole terms we have 



d» 2 (x c ,y c ) 



10 



(21) 



5 kin (x c ,^) = 



3 J dp 



0-2 



d/j d/j' 



[<3(X) - A <3(x c )] A^A^TfAr] 



where d 2 cr 2 / d/j dyu' is defined by Eq. ( IB let and we subtracted the 
term proportional to A = d 2 cr°/ dyu dyu'/( d 2 cr 2 / d/idyu') to cancel 
the temperature-independent contribution, i <2(x c ), caused by scat- 
tering of photons into the line-of-sight. Here d 2 cr!J/ dyu dfi' denotes 
the cross section d _ cr 2 / d/j d/j' for /? = 0. 

In Fig. [6] we show the signals caused in second order of the 
cluster velocity. Clearly, the correction already becomes very small, 
and if relative precision ^ 1(T 4 - 10~ 3 is not needed, these terms can 
certainly be neglected. For the signal related to the motion-induced 
quadrupole temperature corrections are again only visible at high 
frequencies, reaching as 10% for C = 0.05. The contribution to 
the SZ signal from the motion-induced monopole term vanishes in 
the limit 9 C — > 0, since Thomson scattering leaves the monopole 
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Figure 7. Comparison of the different SZ signals for e = 0.05, f} c = 0.01 
and fi c = 1,0,-1 from top to bottom panel. All temperature corrections 
are included for each of the contributions, and the curves were obtained 
using the 3-dimensional integrals for monopole through quadrupole scat- 
tering. We also show the result obtained with SZpack carrying out the 5- 
dimensional collision integral explicitly. Notice that we multiplied the sig- 
nals in second order of /3 C by appropriate factors to make them more visible. 



unaltered. As Fig. [6] also indicates, the convergence of the expan- 
sions, Eq. d 1 4b and dl7t is very slow, and good agreement is only 
reached for low electron temperature. We confirmed numerically 
that the full 5-dimensional collision integral for the monopole and 
quadrupole correction is correctly represented by Eq. (120 ) and ( 12 1 b - 
In Fig. [7] we illustrate the different SZ signals for e = 0.05, 
y3 c = 0.01 and some values of /j c . ForyU c = the kSZ signal vanishes 
and the dominant, but very small, corrections to the thSZ effect are 
second order in f} c . The motion-induced correction to the monopole 
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Figure 8. First few functions Z* m , q™, and S f n for 6» e-0 = 0.03. In SZpack 
we computed up to k = 20. The functions Zj™ are needed to include tem- 
perature corrections related to the motion-induced monopole (oc fi 2 .), while 
C^ ln and Sj™ allow computing temperature corrections to the kSZ (oc /J c // C ) 
and qkSZ effect (oc /J^i^ijUc)). We rescaled by appropriate factors of k and 
fJ c to make them comparable in amplitude. 



is independent of p c , as expected. Also, for p c = -1 the kSZ and 
qkSZ signals add while for p c = 1 they have opposite sign. 



4.1.4 Improved representation of kinematic SZ corrections 



The discussion of Sect. 14.1.31 shows that the expansion in small 
frequency shifts for the kinematic corrections also suffers from the 
same limitation as the pure thSZ effect. This problem again can be 
cured by defining a different set of basis functions to account for 
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Figure 9. Comparison of the expansion given by Eq. 1101 and )23t with the 
full numerical result for reference temperature 9 e fl = 0.03. The agreement 
clearly is very good. We considered /? c = 0.01 and fj c = 1 but found similar 
performance for other values of /i c . 



the various temperature terms: 

1 dp dp' 

N C dVi d k e-^ lff i 



6k\J 

d f dVi 

Ck {X) - Hj dpdp' 



Sf\x) 



N_ f dV 2 
W.J dp dp' 



dp dp' 0*9; 
2 ~ d k e-^ llf i 



dp dp' &*ei 



[M(x') - M(x)] dpdp'n 2 dT], 
g{x')dpdp'rrdr,, (22) 
[Q(x') - A Q(x)] dp dp'rfdrj 



Similar to the thSZ signal (Sect. l3TT~H . we avoided taking the tem- 
perature derivative of the normalization factor, N(0^.) by pulling it 
in front of the integral. Using this basis, the kinematic corrections 
to the SZ signal in the small optical depth limit can be expressed as 



A m;)X 



Ar c g(x c ) + J]Az lk) Cf n (x c ) 



k=0 
i 



(23) 



+ 0i £ Az^Zf (*c) +0 2 C P 2 (J4) -^f0(x c ) + ^ Az (k) Sf\x c ) 

In Fig. [8] we show the first few basis functions for Zj 011 , C]™ 1 , and 
S]™ for illustration. For precise representation of the SZ signal we 
numerically computed the basis up to £ max = 20 with relative error 
set to < 10 -6 . In comparison to the normal asymptotic expansion 
the new set of basis functions again exhibits far fewer oscillations 
even for higher values of k. Equation d23 b therefore allows very 
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AL= pYAf 

Figure 10. Illustration for the transformation of scattering optical depth. The received photons leave the scattering slab of medium at x 1 ^ = *2 and time t^ = tj ■ 
This is the space-time point at which we synchronized the wavefronts. The length of the slab measured inside the cluster frame is L c = x c 2 - x\ = cAf (upper 
panel). A photon that is scattered at Jtj at time rj = - Af into the line-of-sight is received by the observer at the same time as a photon that is scattered at 
£ and time £. — When the slab is moving with velocity v = cfi towards the observer (lower panel), the momentary length of the slab is reduced to L = L c /y 
by Lorentz-contraction. A photon that is scattered at jq(fi) = X2te) ~ + ^0 = cAt is received by the observer at the same time as a photon scattered at 
X2(fl) = jfClf) towards the distant observer. During the interval At the end of the slab moved by AL = ficAt to X\i)% = t^) = X2 — L, We therefore find the 
effective time-interval At = L/c(l — f£) = Af /y(l — 0) = y(l + /}) Af° over which the scatterings of photons that are received simultaneously by the resting 
observer take place. For illustration we chose /? = 0.6, so that L = 0.8L C , AL = l.2L c and At = 2Af. 



precise computation of the SZ signals for an observer at rest in the 
cluster frame. 

With Eq. J23b we can now compute the SZ signals in a very 
fast and economic way. In Fig. [9] we illustrate the performance of 
this expansion. For all considered cases we find agreement at the 
level < 0.001% with the full numerical result in second order of y8 c , 
at least away from the crossover frequency, where the total signal 
is small. For temperatures # e < 0.01 one can use the asymptotic 
expansions given above to compute the SZ signals. A combination 
of both set of basis functions therefore allows calculation of the SZ 
signal with very high precision. We implemented these approxima- 
tions in SZpack, and found that even for # c as 0.06 - 0.07 the new 
set of basis functions still represents the full numerical result with 
precision < 0. 1%. We also checked the expressions for other values 
of fi c and found excellent agreement with the full numerical result. 

To emphasize the improvement in the performance, evaluation 
of the full 5-dimensional collision term on a standard laptop takes 
about 2 minutes for 40 frequency points with SZpack, while calcu- 
lation of the SZ signal with our improved set of basis functions is 
practically instantaneous. For computations of the SZ signals from 
simulated clusters this makes an important difference. 

4.1.5 Additional small corrections 

When considering terms in second order of yS c , in reality one should 
also include the primordial CMB quadrupole, octupole and other 
large-scale anisotropies, which all have a spectrum x c d y jip\(x c ). For 
clusters with /? c a 3 x 10~ 3 these in principle have comparable am- 
plitude (^ 10~ 5 - 10~ 4 ), and in different orders of the electron tem- 
perature the anisotropies they do couple to the Compton scattering 
cross section. With the 5-dimensional collision term we can esti- 
mate this effect using SZpack. We found that the correction always 
entered at the relative level of An/n < 10~ 4 , except for close to the 
crossover frequency, where the correction could reach a 0.1% of 
the kSZ signal. However, in absolute terms this correction is tiny. 
Similarly, higher order motion-induced corrections remain small. If 
extremely high precision was required it would be straightforward 
to include these terms. We also confirmed that the effect of stimu- 
lated scattering and electron recoil become important at a relative 
level An/n < 10~ 6 . 



4.2 Transforming the SZ signal into the CMB rest frame 

The expressions for the SZ signal discussed so far are only valid 
for an observer at rest in the cluster frame. Also, they describe 
the change of the photon occupation number caused by electron 
scattering assuming the optically thin limit, and multiplication 
by 2hv 3 /c 2 gives the distortion of the spectral intensity signal. 
As mentioned above, in the cluster frame the derivation is rather 
straightforward and very clear in terms of the required scattering 
physics. However, the transformation to a general observer be- 
comes slightly more complicated, especially when thinking about 
real measurements and the appearance of the cluster to the observer 
(see Sect.l4~5l 

Below we consider the transformation of the SZ intensity sig- 
nal into an arbitrary observer frame. However, to understand the 
problem it is easiest to start with the special case of an observer 
at rest with respect to the CMB, where the CMB spectrum is given 
by a blackbody, ;ipi(x), with no anisotropies. From the cosmological 
point of view this is also the most relevant frame, for instance when 
discussing about the statistics of the large-scale velocity field. The 
first obvious transformation relates the frequencies of the observed 
photons: .y c = y c x{i — yS c /j c ), accounting for the Doppler effect. 
Transformation of y c gives nl = (ji c -j3 c )/(l -/? c /u c ), describing 
the relativistic aberration effect. Here ji c = f) c ■ y is measured in 
the CMB frame, while the azimuthal angle under which the scat- 
tered photon is received does not change. These expressions can be 
directly inserted into the formulae describing the SZ signals, e.g., 
Eq. J 1 0b and l !23t . when using our improved set of basis functions. 
However, there are two additional transformations that deserve a 
few comments. The first is related to the optical depth and the other 
to the electron and photon temperature. 

4.2.1 Transformation of the optical depth 

Let us start by considering the optical depth along a given line- 
of-sight. Because time along the p hoton's path transforms aj* 2 ! 
(compare ISazonov & Sunvaevl |l998) dr c = y~\\ +/3 c (4)- 1 dt = 
Tc(l - j8 c ^c)df, for the differential of the optical depth we have 

12 We denote time with t and spatial coordinate with r in the CMB frame. 
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&T°(f,T*) = y c (l - ct i c )ccr J N c t lf,r e )<it = (1 - fcn c )dr(t,r) 
where we used N^if, r c ) = N e (t, r)/y c (also see IChluba et alj2012L 
for more detailed explanation). The Lorentz factors cancel out af- 
ter combining time-dilation and length-contraction. The remaining 
factor of (1 - yS c yu c ) accounts for the effect of retardation: while 
photons are propagating along a given line-of-sight the cluster is 
moving with f3 c p c along that same line-of-sight. For a fixed inter- 
val At = Af this shortens the apparent distance that photons travel 
through the cluster medium by (1 - f3 c p c ). 

However, in terms of the SZ signals, the factor (1 - y0 cy u c ) 
should not be considered separately, since physically the rele- 
vant quantity is the total number of scattering events, iV sc (y c ), 
along the photon path through the ICM, which is Lorentz-invariant 
and related to the integrated scattering optical de pth (see e.g., 
iRvbicki & Lightmanlll979l ; lAbramowicz et aill99ll) . Different ob- 
servers therefore always agree on the number of scatterings encoun- 
tered by the photons along a given geodesic through the cluster. 
The apparent direction under which the photon is received along a 
particular geodesic and also the specific moments of the scattering 
events depends on the frame because of the aberration effect, time- 
dilation and also geometric effects, but Af sc (y c ) is always identical. 

To understand this point better, let us consider a simpler case 
in which a thin slab of scattering medium is moving with yS along 
the jc-axis towards the observer (see Fig.[H)]for illustration). If the 
cluster and observer are at rest with respect to each other, then we 
have the optical depth At c = J^' 2 ccr x A^ dt c = cr T A^L c cc N sc . Here 
we assumed a constant number density, N%, of electrons along the 
slab. As indicated, At c is directly proportional to the number of 
scatterings along the photon path, or equivalently, the number of 
electrons encountered by the photon along the given world line. 

Using the transformation for dr given above, we also have 
Ar £ = (1-/3) f 2 ccr T N,,dt = (1 - yS)ccr T A' c Af. As illustrated in 
Fig. [TO] the effective interval over which the scatterings take place 
increases to At = Af ly{\ - (3) = L c /cy(l -/?). We therefore have 
(1 - j3)ccr T N e At = o-jNelS/y = At c , because iV£ = N e /y. More 
generally, this implies that the integral Ar*(y) = f t ~ ccrjN e dt has 
the interpretation Ar'(y) = At c (7 c )/(1 — yS c /j c ) in terms of the clus- 
ter frame line-of-sight optical depth along the photon's world line 
defined by y c . This shows that Af(y) is no longer just related to the 
number of scatterings (or the scattering optical depth in the classi- 
cal sense) along the photon's path but includes a kinematic factor 
(1 - ficHc) 1 to account for the change of the electron flux in the 
CMB rest frame. This implies that the variable At* (y) describes a 
mixture of scattering effects and kinematic correction terms. 

To compute the SZ signals in different observer frames we 
directly use AT c (y c ), only transforming the apparent directions un- 
der which the photon are received. This variable is very convenient 
because At c oc N x without any extra dependence on the cluster's 
(or observer's) motion, such that a clean separation of scattering 
and kinematic effects is achieved. As we explain in Sect. 1 4.41 in 
comparison with previous works (e.g.. ISazonov & Sunvaevlll998l ; 
iNozawa et allll998ll200^ : lshimon & Rephaelill2004l) this makes a 
difference for the interpretation of the kinematic correction terms 
at order 0(J3 c 9 e ) and 0(fl 2 c ). 

4.2.2 Transformations of electron and photon temperature 

The effects mentioned in this section only become relevant for fu- 
ture high, precision SZ observations, because they enter the SZ sig- 
nal at order 0{fi c B c ) and O(0l0 c ). Nevertheless, they touch on inter- 
esting aspects connected with the determination of the cluster frame 



electron temperature, 7£, and the CMB temperature, T r , which are 
worth discussing briefly. In particular, when considering the prob- 
lem of eliminating from the expressions for the SZ signal, e.g., 
with X-ray observations, this section should be of interest. 

It is evident that only in the rest frame of the moving cluster 
(or a small volume element with constant temperature and number 
densi t y) it is possible to define 7^ in a meaningful way (see [Reiser] 
1 1994 Ichacon-Acosta et alj|2010l , for some related discussion). In 
this sense the temperature of the gas becomes a parameter that, 
like the rest mass, describes the properties of a parcel of gas. An 
observer, who is moving with respect to the cluster frame, can only 
indirectly infer the rest frame electron temperature, for example, by 
measuring the widths of thermally broadened resonance lines, the 
shape of the X-ray spectrum, or even by measuring the SZ signals 
discussed here. In all these cases we observe some photons and 
by using the Lorentz-invariance of the photon occupation number, 
after making assumptions about the physical processes causing the 
signal, we can in principle understand how to infer the electron 
temperature defined in the rest frame of the cluster. Consequently, 
the 'measured' value for the electron temperature, T' t , is directly 
linked to the method that was used to obtain it. 

Thinking, for example, of a thermally broadened resonance, it 
is clear that to first order in yS c only the position of the line changes, 
but its fractional width remains the same, so that one can directly 
infer r c ' « T£. In second order of the motion this is no longer true, 
and the detailed transformation law for the photon distribution into 
the observer frame has to be considered. When using the cluster's 
X-ray spectrum to determine the electron temperature, it is also im- 
portant that to first order in f3 c the Doppler effect again changes the 
received photon frequency. This implies that the inferred value of 
the electron temperature (by fitting the shape of the X-ray spec- 
trum) is T' t =: TH (I - /JcA'c)- However, in both the narrow line and 
X-ray spectrum case degeneracies with the redshift of the cluster 
exist. This shows, that as long as kinematic corrections are insignif- 
icant (say in comparison to the measurement errors) we need not 
worry about the way T£ was inferred; however, small differences do 
arise at high precision, and dependencies on the method becomes 
an issue. Observationally, measurements of the electron tempera- 
ture (e.g., from X-rays) at the sub-percent level are futuristic, so 
that this effect will not become significant any time soon. 

For similar reasons, determinations of the monopole temper- 
ature of the CMB, T y , must be considered with care. In the CMB 
rest frame, by averaging the spectrum over the whole sky, we di- 
rectly find the value of T y , by comparing the average spectrum 
with a blackbody. The same result is of course obtained by mea- 
suring the spectrum in any single direction (we neglect primordial 
CMB anisotropies). However, for an observer who is moving with 
respect to the CMB, the 'observed' value for the CMB temperature, 
T' in second order of the peculiar motion again depends on the 
observational strategy. By subtracting the motion-induced dipole 
and quadrupole temperature anisotropy in the moving frame, the 
mono pole temperature is T' ss 7" r (l -0^/6) dChluba & Sunvaevl 
120041) . If on the other hand, we just integrate the observed inten- 
sity (including the motion-induced dipole and quadrupole spectral 
anisotropies) over the whole sky, we find an average energy density 
of p y x a^Ty{\ + |y8j) because of the s uperposition of blackbod- 
ies with spatially v arying temperatures dChluba & Sunvaevll2004l : 
IChluba etal]2012l) . Here a R is the radiation constant. We therefore 
can infer an effective temperature of T' as T y (l + 01/3). 

These are just two simple (though not necessarily practical) 
examples of how the measurement procedure determines the rela- 
tion between the observable, T' y , and the rest frame CMB temper- 
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Figure 11. Kinematic corrections from a change of the observer's frame. 
The solid/black line shows the SZ signal for a cluster with fl c = 0.03, 
y3 c = 0.01 and fi c = 1 in the cluster frame. The other two lines show the 
correction to this signal (times 10) when transforming into the CMB rest 
frame (dashed/blue line) and for an observer who is moving in the opposite 
direction of the cluster with j8 D = 0.01. 



ature, T y . Since in our expressions derived for the SZ signals both 
Tl and T y appear, it is important to bear the above comments in 
mind, when interpreting measurements at high precision. The val- 
ues of T and T' e obtained by some method in the moving frame 
affect the inference made for the SZ effect. For T y , this only enters 
the problem via terms connected with the thSZ effect, e.g., Eq. tlOt , 
since all kinematic correction terms are already at least of first or- 
der in j8 c . However, for Tl kinematic terms can even give rise to 
corrections oc B c 8 e . In other words, errors on the inferred values for 
T y and 7£ will become important when interpreting the SZ signals 
with high precision. For the current precision in the value of the 
CMB monopole, this limits the interpretation of SZ signals to no 
better than 0.1%, but in the future the monopole of the CMB 
might be measured with m uch higher accuracy, for example using 
PIXIE dKogutetalJl201lh . 



4.2.3 Final SZ signal in the CMB rest frame 

With the comments above it is now straightforward to obtain the 
appropriate expressions for the SZ signal in the CMB rest frame. 
Using Eq. J X Oi l and Eq. d23l l. for a given observing frequency, x, 
electron temperature, 7£, cluster velocity, B c and direction, p c , and 
line-of-sight optical depth, At c , we must evaluate the expressions 
at frequency x c = y c x(l -B c p c ) and p? Q = (p c —B c )/{1 — B c p z ). In this 
way we obtain a precise prediction for the associated SZ signal up 
to second order in B c . For simplicity one could insert pi a p c and 
x c sa x for all terms that are already of second order in B c , but this is 
not required, as the 'error' enters at 00 s .). Similarly, for terms that 
are already first order in B c one can use pi sa p c - /J c (l - pi) and 
x c ~ x(l -B c p c ). 

In Fig.[TT]we illustrate how the transformation back into the 
CMB rest frame affects the observed signal. The main effect is a 
shift of the photon spectrum towards higher frequencies. For the 
considered case, the correction is comparable to a 5% of the SZ 
signal in the cluster frame. In particular, close to the crossover fre- 
quency, x cr =! 3.83, the kinematic correction makes a large differ- 
ence, shifting x ci towards higher frequencies. 

When looking at an extended object, there is a collection of 
sight lines through the object. Each of them has a different value of 



p' c . Fixing the axis toward the object and assuming that it is mov- 
ing as a whole in one direction with velocity vector defined by p c 
and C relative to this axis, p' c is given by p' c = p c p + cos(0 c - 
<P) vO ~ A<c)(l _ l* 2 )- Here p = cos(0) and <f> defines the line-of- 
sight relative to the reference axis. Typical clusters have angular 
size of only a few arcminutes few x 10~ 4 rad), so we can ap- 
proximate p! Q a i_i c (l - ±g 2 ) + cos(0 c - ~ p 2 8. At a relative 
precision of a 10~ 4 - 10 3 one can therefore neglect higher order 
terms and just use p' c « p c for each line-of-sight through the clus- 
ter. This is equivalent to the flat-sky approximation and the optical 
depth in this case is just a function of position on the sky At c (8, <f>) 
with roughly one common value of p c independent of <f> c . 



4.3 Transforming into the frame of a general observer 

Following the discussion of the previous section, it is now trivial 
to give the transformation laws for a general observer. If the ob- 
server is moving at a velocity J3 a with respect to the CMB, we 
can first transform the SZ signal from the cluster frame into the 
CMB frame and then subsequently relate the observed photon fre- 
quency x„ and photon direction y to the corresponding quantities 
in the CMB frame. To study the large-scale velocity fields we are 
interested in the values /? c and p c inside the CMB frame, and it is 
therefore unnecessary to express these 'parameters' in the observer 
frame. Thus, we can simply use the transformations 



x c = y c x[l -B c p c ] = y c y x [l - B c p c ][l + B p a ], 
Pc-/3 C Po+Po 



l-'c 



1 -B c p c 



H ■ 



1 + B p a ' 



(24) 



with the Lorentz-factors y c = 1 / -^1 - B\ and y = 1 / -Jl — and 
direction cosines p c = j8 c • y, p = fi a ■ y, and p a = j8 Q • y . The 
photon frequency, x , and the value of p are directly measured 
in the observer's frame, but the dependences of the SZ signal on 
the cluster's and observer's peculiar motions and the relevant di- 
rection cosines are still expressed in the rest frame of the CMB. As 
mentioned above, this is a reasonable choice, since in this way the 
interpretation of the measured SZ signal (in the observer's frame) 
directly reveals B c and p c in the CMB rest frame, the quantities we 
are after from the cosmological point of view. 

The value of B can be directly determined by looking at the 
CMB temperature dipole inside the observer's frame. For our mo- 
tion with respect to the CMB rest frame we have B = 1.241 x 
1(T 3 (1 ± 0.2%) toward s (l,b) = (264.14° ± 0.15°, 48.26° ± 0.15° ) 
in galactic coordinates dFixsen et alJl996tlFixsen & Mathej2"002T) . 
The value of p a is then also fixed, by comparing the position 
of the cluster on the sky with the dipole axis. Since B is small 
the additional modification caused by our motion relative to the 
CMB re st frame also r emain s small. This effect was first investi- 
gated bv lChluba et"ai] < l2005h . T heir expressions for the SZ inten- 
sity change were confirmed by ( INozawa et al.ll2005l) . It was also 
pointed out that the motion of the Solar System s hould lead to 
a dip olar anisotropy in the cluster number counts dChluba et all 
120051) . Furthermore, the small errorbars on B a and the direction of 
the dipole imply that the effect of our motion with respect to the 
CMB can be taken out with ^0.1% precision. 

With Eq. I l24t . we can verify that x c = x for j8 c = j8 . Here it is 
important to note that in the observer's frame we measure p = p c c , 
and hence (1 - B c p c )(\ + B p B ) = (1 - B c p c )(\ + B c p c c ) = yf. For 
an observer at rest in the CMB frame we also find x = x and 
p = p, as expected. Note that by using the expressions i24\ we 
can avoid additional approximations related to the transformation 
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into the observer's frame. The precision of the predicted SZ signal 
is therefore only limited by how well we are able to represent the 
distortion in the cluster frame. 

In Fig. Qj] we illustrate how the SZ signal is affected by the 
transformation to an observer moving with a rather large velocity 
in the opposite direction of the cluster. For comparison, our motion 
relative to the CMB rest frame causes a maximal correction that is 
about one order of magnitude smaller. Notice that for the computa- 
tion of the SZ signal we used our expansions, Eq. ( 110) and Eq. d23l >, 
and only needed to fix At, T e , T y , f3 c , p c , /?„, and p. . Here, we as- 
sumed that the cluster is isothermal along the line-of-sight and we 
omitted bulk flows of the ICM. From CMB observations we can di- 
rectly infer T y ,ji , and p . The remaining parameters describing the 
SZ signal are At, T c ,f3 c , and /j c ; these can be directly constrained by 
multi-frequency SZ observations, where At and T e now denote the 
cluster frame scattering optical depth and electron temperature, re- 
spectively. Generalization to the non-isothermal case and inclusion 
of bulk flows of the ICM is straightforward and will be discussed 
in our future work. 



not explicitly addressed previously, but with t he explanations given 
in Sec t. l4.2.1 | it is straightforward to sh ow that lSazonov & Sunvaevl 
dl998l) and lShimon & Rephaelil J2004l) expressed the SZ sign al in 
terms of At* =_At/(1 - /3 c p c ), while iNozawa et al.l i 19981) and 
iNozawa et al .1 l [2006h used At" = At* /y c . These redefinitions ex- 
plain the above difference and consequently imply alternative in- 
terpretations of the kinematic corrections to the SZ signal in each 
case. From the physical point of view the variables At* and At** 
both describe a mixture of scattering and kinematic terms, while 
with our choice of variables these effects are cleanly separated. 

Similarly, for the kinematic term oc 1 (with no additional de- 
pendence on yu c ) we find 



A/f 



llksz (;r,r) 1 , „ 



3ff(x)] = -p 2 [Y (x) + g(x)] (27) 



using D* x @ = -Q(x) + 3G(x) = -M(x). This dif fe rs by -§G(x) 
from the results of ISazonov & Sunvaevl dl998l>. I Nozawa et al.l 
dl998l) . IShimon & Rephaelil J2004I) . and INozawa et alj J2006I) . for 
the same reason as above. 



4.4 Comparison with previous results (CMB rest frame) 

We now compare our results for the kinematic corrections directly 
with the expressions given by different groups. For this we consider 
the SZ signal as measured in the CMB rest frame. Following the 
discussion of Sect. 14.21 and using the analytic expressions for the 
asymptotic expansion given above, for f3 c <K 1 we find 



An(x, y) 

At(t>) 



Yf + \o x Y k 
b 



g + e^^l^-xdj,} 



(25) 



We emphasize again that At and r e denote the cluster frame scat- 
tering optical depth and electron temperature, respectively. All 
frequency-dependent functions are evaluated at frequency x and 
we used p? = f/^O") + *• We also defined the differential oper- 
ators T> x = 3xd x + x 2 d\, D' x =2 + xd x , and £>** = 1 - xd x , and 



{& + ft Yj k k Zt %Dk B )- As before, we as- 



the spectral function "K = l@ + e 2J; n ™ ' 

sumed that the cluster is isothermal and has no internal bulk flows. 

With Eq. l |25t we can easily read off terms of different orders 
in the electron temperature and f} c . As mentioned earlier, our ex- 
pression for the purely thSZ effect agree wi t h those given earlier 
in the literature fchallinor & Lasenbvl 1998; Sa zonov & Sunvaevl 



ll998l : lltohetanil998l ; IShimon& Rephaelil |2004|) . Also the kSZ 
term, oc fi Q p z Q, is the same. For the qkSZ effect we however find 



AT(r) x 30&^« 



(26) 



in the CMB frame, where we used T)**Q = Q(x) = Q(x)X(x). 
This result differs by — ^0jP2(Uc)&( x ) fr o m the correspondin g 



:orrespoi 

terms given ir ^1 Sazonov & Sunvaev ( 1998), Nozawa et al.l JT998T) 



IShimon & Rephaelil ( |2004|) . and lNozawa et alj J2006I) . 

The difference arises because other conventions for the opti- 
cal depth variable were used in these earlier works. This aspect was 



13 IChallinor & Lase nbv 1 1999) only gave the kinetic equation to describe 
the scattering of CMB photons by moving electrons, but did not integrate it. 



4.4.1 Temperature and velocity cross-term in first order of y8 c 
From Eq. &25) in first order of the electron temperature we obtain 



An x (x,f) 
At(7) 



Beetle [D^-xd x Y Q ] 



9^ 



138 42 7 
-— &-— Y Q + -xQm + x) 



(28) 



Again this expression differs from the o ne of ISazonov & Sunvaevl 
ll 19981) and INozawa et all dl998l l2006h by the amount expected 
from the redefinitions of the opti cal depth variables explai ned 
above, i.e., -O^p^fi^Yo. The result o f Sh imon & Rephaelil J2004 for 
this term differs from ours by additional -Qji^n^Q. We attribute this 
to a slight inconsistency related to the use of th e Lorentz-invariance 
of the electron phase space distribution in Sh imon & Rephaelil 
d2004l) , but did not investiga te this question any f urther. This dif- 
ference was also discussed in INozawa et alj feoOrj) . 

To conclude, at first order in /3 C we find full agreement of 
our results for the asymptotic expansion of the Boltzmann colli- 
sion term with those given in pre vious works dSazonov & Sunvaevl 
ll998MNozawa et alj 19 98. 2006); however, only after including the 
extra factor of (1 — /? c/ u c ) caused by the retardation effect in the defi- 
nitions of the scattering optical depth. This factor is absorbed by the 
line-of-sight integral and does not appear separately, as explained in 
Sect. |4.2.T| When interpreting the SZ signals of individual clusters 
this difference can become important at the level of a few percent. 



4.4.2 Temperature and velocity cross-term in second order of{3 c 

Fo r similar reasons, the tempe ratu re correction oc 2 d s P2(Lt c ) given 
bv lNozawa et~aT] d 1 99 8L [2000) and lShimon & Rephaelil J2004l) . and 
all other second order kinematic terms also differ from our expres- 
sions by some small amount. For example, our first order tempera- 
ture correction oc f} 2 6 c P2(p c ) reads 



An^ sz - m (x,f) 



Ar(r) 



~0te s p 2 (Mo) 

128 



Q^ + l(\- x d x )D^ + ^x 2 d 2 x Y a 



g+^-Y + 6x 3 d x ii PI + ^x 4 d x n PI 



(29) 
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where we used Eq. ( 1 1 8t to simplify the expression. The additional 
term related to the retardation factor (1 - /? c p c ) from the optical 



depth gives rise to -(2/3)/?^ 
Ani ksz - (1) 0,r) 



At* (7) 



■PiiUcXD^-xd.Yo), yielding 
44£ + 12y + 6x 3 d> P1 + ^^"pi 



With E q. dAlt one can show the equivalence to the expression of 
iNozawa et al.l dl998h . noting that At* = At" + 001). However, in 
the present work the factor (1 - /? c p c ) is not considered separately, 
but as part of the optical depth integral. Otherwise, the interpreta- 
tion of optical depth is different and does not only reflect the effect 
of s cattering, as explained above. 

INozawa et alj fl998[) did not explicitly carry out the transfor- 
mation of the electron number density into the CMB rest frame. 
This introduces another factor of y" 1 ~ 1 - jfil, which affects the 
expressions for the terms oc 0^6^ (without extra /^.-dependence). 
For example, we find 



A „mkSZ,(l)( Xj ^ 



At( 7 ) 



2Y^--(2 + xd x )^ 



234 r+ 66 v x 77 W x 7 4.4 

— y + y TE X x 15 x x 



where we made use of Eq. OS), and D x Y k = 6Yf n . Transforming 
to At* gives an extra term, —^0^d e W^ n — xd x Yo], so that 



An mksz - (l) (jf, y) 
At*( 7 ) 



1 2 

3*0 + ^ 



This is prec i sely \pl6 c Y(, smaller than the expression of 
INozawa et al. I Jl998h ." which, as mentioned above, is because N% 
was used to define the optical depth integral instead of N e = N%/y c . 
However, this modification is small and should not affect the final 
result and interpretation of the SZ signal significantly. We also note 
that the difference in the definition of N e does not affect the lowest 
order ji\ term, Eq. i BOI l. since Thomson scattering leaves the CMB 
monopole spectrum unaltered. 



4.4.3 Energy-transfer and conservation of photon number 

From the final solution for the cluster SZ signal one can compute 
the change in the number density of photons and also the total 
Compton energy transfer caused by th e moving cluster. This can 
directly serve a s a consistency check dSazonov & Sunvaevlll998l : 
lltoh et alii 19981) , because the total number of photons should not 
change and also the energy transfer rate is well-known in this case. 
However, when considering different lines-of-sight through a ex- 
tended object and the average properties of the received photon in- 
tensity, some subtle effects arise which we explain here. 

Before carrying out the line-of-sight integral, from the photon 
Boltzmann equation we have 



dn(x,y) , An(*.f) 

~ (1 - j3 c p c )c(T T N e 



dt 



AT(y) 



(3D 



where An(x ,y)/ Ar(y) is given by Eq. l|25l >. This expression is con- 
sistent with lSazonov & Sunvaevl l ll998l) for all orders in fj c and 8 e 
presented there. Up to second order in /3 C , it is also equivalent to the 
expression obtained fro m the Lorentz-boosted Boltzmann collision 
term, which was used bv lNozawa et al. 1 1 998[) an d more recently in 
the cosmological context by Chluba et al. 1 20121). However, ^ex- 
plained above the collision term presented bv lNozawa et al. i mi) 
has to be multiplied by l/y c in addition, to transform N% to A^. 
For a fixed time, from Eq. ( 13 U we can carry out the angle 



averages and compute how much the photon field is affected by 
scatterings at a given location. Since \ j Pi(M)dM = for > 0, 
and because \ j pP\ (p) dp = | and ^ j pPiip) dp = 0, we find 

(^)"*£<(* + W)-S*»S"<W 



Considering, for example, the case 



k=0 

we have 



xd x Y k \ 



I dn(x,y) 
\ ccr T N c dt 



n L .=() 



a l -PlQ(x)-^plg= l -plY (x). 



(32) 



(33) 



This expression implies ( — p ) 



for the change of the local 



* /e c =o 

photon number density, confirming that the scattering event con- 
serves photon number. We also find ^ ^ j _ Q = |y8^p r ccr T A' c 
for the change of the energy density, which is in agreement with 
the classical result for the Compton energy transfer (e.g., see 
iBlumenthal & Gould [l97ol : iRvbicki & Lightmanlll979l) . We con- 
firmed that higher order terms also yield consistent results with re- 
spect to photon number conservation and energy transfer. 

However, an important difference arises if we first carry out 
the line-of-sight integral and then average over different photon di- 
rections. Let us again consider the terms for 6> e = 0. From Eq. l !25t 
and d27t we have 



A«(x, y) 
At( 7 ) 



&9, 



which means 
I AA^ 
~At 



and 



8c =0 



At 



« L .=o 



X 3# P - 



(34) 



(35) 



Not only does this apparently violate the photon number conser- 
vation, but also the transferred power is two times larger than that 
given by the classical formula. What is happening here? 

The explanation is related to the differ ence between emitted 
and r eceived power from a moving source dRvbicki & Lightmanl 
[1971) . For the classical expression given above, the power emitted 
at a single point is integrated over all directions for any given in- 
stance of time. Imagining several sources embedded along a slab of 
medium, the average emitted power from all these sources is pro- 
portional to the number of sources inside the slab. However, the to- 
tal power emitted by these sources was computed for each individ- 
ual source separately and then averaged over the different sources. 
This is equivalent to placing all the sources into a single point and 
then averaging the emitted photon field using time-intervals mea- 
sured in the source frame. In contrast, the averages Eq. (1351 are 
carried out such that the time-intervals are measured by the ob- 
server, with the retardation effects consistently included. For SZ 
observations by some distant observers placed on a sphere around 
the cluster, Eq. l !35t is therefore the relevant expression. 



4.5 Appearance of the SZ cluster and additional kinematic 
effects on the measured quantities 

Thus far we have only considered the transformation of the photon 
occupation number and how the spectral intensity along a given 
line-of-sight is affected by kinematic corrections. However, addi- 
tional effects arise in practice. For example, due to the aberration 
effect different rays of photons change their separation: in the di- 
rection of the motion they become closer, while in the opposite 
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direction they diverge from each other. This means that when ob- 
serving the flux of photons in a given direction and some fixed solid 
angle in different frames, the signal changes accordingly. Conse- 
quently, the apparent angula r size of the cluster changes by a small 
amount JChluba et alj2005h . depending on the direction of the mo- 
tion of the observer relative to the cluster. Furthermore, the geom- 
etry of the cluster affects the apparent shape on the sky, and devi- 
ations from non-isothermality and the radial distributions of elec- 
trons modify the SZ signals. If we were to compare the results ob- 
tained by different observers or in different directions on the sky, 
this effect must be taken into account. 

Also, in observations the photon flux is integrated over some 
frequency filter and with a finite sampling rate. Fixing these in dif- 
ferent frames implies kinematic corrections to the 'observed' SZ 
signals for different observers. Similarly, it matters if photon flux 
or intensity are consider. We will discuss these aspects of the prob- 
lem in more detail in another work. 



4.6 Astrophysical contaminations 

There is a long list of additional effects that affect the SZ sig- 
nal or compromise its interpretation at some significant level. For 
example, we neglected the effect of multiple scattering (see e 



example, we neglected the effect of multiple scattering (s ee e.g . 
Rephaelilll9953; iMolnar & Birkinshawlll999l; [DoTgov et alj|200ll ; 



Shimon & Rephaeli 2004 : Nozawa et al J 12003) . which should be 



most noticeable close to the crossover frequency. It is also known 
that the ICM should host populations of non-thermal electrons 
(e.g. lEnfilin & KaiseJl200d ; lBlasi et~ai]|200Cl : IColafrancesco et alj 
120031) . These give rise to non-thermal SZ distortions, which we 
omitted here, but do affect the SZ signal, especially at high frequen- 
cies. Also, the local radiation field inside the cluster is not neces- 
sarily described by a pure blackbody but could be contaminated, 
e.g., by synchrotron radiation. In that case one must also consider 
the effect of up-scattering for these parts of the photon distribution, 
another problem we leave for some future investigation. Moreover, 
galaxy clusters can act as a moving lens, deflecting photon from the 
CMB background and causing a characteristic temperature pattern 
whic h can reach a level comparable to ^ 10% of the kSZ e ffect 
(e.g.. iBirkinshaw & Gulllll983l : iMolnar & Birkinshawll2003l) . Fi- 
nally, we did not account for the effects of the inhomogeneous ICM 
structures and a su ite of astrophysica l processes, such as turbulent 
gas motions (e.g.. [Dolag et alJl2(5o3 : lLau et al J [20091: IVazza et alj 
2009MNelson et alj201 11) or magnetic fields (e.g.. lKoch et alj2003l ; 



Gopal & Rovchow dhurv 201 3), on the SZ signal. For accurate pre- 



dictions and inteipretation of the SZ signals all these processes 
must be included. 



5 SUMMARY 

We considered the SZ signals caused by the scattering of CMB 
photons by hot electron residing inside clusters of galaxies, ob- 
taining improved expressions that allow fast and precise computa- 
tion of these signals. We compute the effect in the cluster frame 
and then derive the signal for the general observer by Lorentz- 
transformation. This procedure is physically well-defined and al- 
lows us to understand the origin of the different kinematic correc- 
tion terms in a simple way. The method provides a quasi-exact 
representation of the Compton collision term, with frequency- 
dependent basis functions that are informed by the underlying 
physics of the scattering process. The precision of the calculation 
is limited by the accuracy achieved inside the cluster frame, which 



we push well beyond the precision of current SZ observations and 
uncertainties associated with cluster physics (see Sect. 14.6b . 

As explained here, the kinematic corrections to the SZ in- 
tensity signal obtained in previous analyses differ from ours by a 
small (comparable to a few percent for a typical cluster) amount 
(Sect. |4~4t . This is because here we explicitly express the SZ sig- 
nal in terms of the cluster frame optical depth along the photon's 
path. This choice of variable provides a clean separation of kine- 
matic effect from those related to scattering events. The difference 
might be important for future precision measurements of the clus- 
ters bulk velocity and when confirming the redshift-scaling of the 
CMB temperature using the SZ effect. In Sect. [43] we also briefly 
mention additional kinematic effects on the appearance of the SZ 
cluster and related instrumental issues; however, a more detailed 
analysis is beyond the scope of this paper. 

One of the main product of this work is SZpack, a numeri- 
cal library which allows fast and precise (< 0.001% at frequencies 
hv < 20kTy) computation of the SZ signals up to high electron tem- 
perature (kT c 25keV) and large peculiar velocity (v/c - 0.01). 
SZpack should be useful when analyzing the SZ signals of high res- 
olution, multi-frequency observation carried out, for exam ple, with 
ALMA dTarenghilhooljh and CCAT dSebring et al]|2006l) . Further- 
more, the package should enable fast and accurate calculations of 
the predicted SZ signals from high-resolution clusters simulations. 
Extensions to account for the effect of non-thermal electrons us- 
ing the method described here are planned. Similarly, ou r approach 
should be applicable when p olarization effects (see, e.g. Jltoh et alj 
120001 : 1 Challinor et alj [20001) are considered. The impact of other 
processes mentioned in Sect. 14.61 will be studied more carefully in 
our future work. 
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APPENDIX A: COMPUTING DERIVATIVES OF THE 
PLANCK SPECTRUM 

The derivatives of a Planckian distribution can be represented in 
different ways. To compute the temperature corrections to the thSZ 
effect up to /7 th order of 6 e we need derivatives up to ordei0 



14 To analytically include all kinematic corrections in O(fil) two additional 
derivatives are needed. 
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k = 2n+2. Here we give explicit expressions for k < 8; however, we 
also show how to compute these derivatives up to any order using 
Eulerian numbers. This procedure is useful, since the expressions 
are rather lengthy while the relations using Eulerian numbers are 
basically analytic, involving summation of previous coefficients. 
Also, the recursions given here greatly increase the numerical sta- 
bility for the evaluation of higher order derivatives. 

Al Explicit expressions up to the 8 th derivative of n Pl 

In lltoh et alJ Jl998h and lSazonov & Sunvaevl dl998h the derivatives 
jc*(9*/ipi were expressed using the variables 



X( X ) : 



*coth(i), 



S (x) = xj sinh 



(!)■ 



To find the corresponding expressions for different derivatives of 
the Planck function, «pi(x) = 1 j{e x - 1), the relations 

x k+l d k l +l = (xd x -k)x k d k x 
xd x X = X- -S 2 

xd x S = S - -XS 

are very useful. For example, in terms of X and 5 one can write 

xd x nn = x 2 d 2 x n Pl =@X, 



x 3 <9 3 n P , = -Q 



x 2 + i§- 



x*tf x n PX =gx \x 2 + 2S 2 ] 



where in addition we introduced the function 

xe x 



§{x) 



S 2 



, v n, - xn n (x)\\ + n Pl (x)] = —. 
(e x - 1) 4x 



(Al) 



(A2) 



However, the above expressions by no means are unique. For ex- 
ample, one can also write x 3 3 3 ra P i = -£?^5 2 + X 2 ]. Since Q is 
directly related to 5 we therefore tried expressing everything with 
just X and Q. Using the relations 



xd x Q = Q- QX 
xd x X = X- 2x@ 
X 2 = Ax@ + x 1 
it is straightforward to show that 
xd x n P1 = -Q, 
x 2 d 2 x n PI = +gx 
x 3 d 3 x n Pl = -xQ [6@ + x] 
x 4 dtn Pl = +xQX[\2g + x\ 



(A3a) 
(A3b) 
(A3c) 



x 5 d x ii Pl = -xQ \30QX 2 + x 3 ] 

x 6 tf x ti Pl = +x 2 gx [m@ (<sg + x) + x 2 ] 

x^lnpi = -x'Q [A2Q (40Q [3Q + x] + 3.x 2 ) + x 3 ] 
x s dln Pl = +x 2 QX [252ff (20ffX 2 + x 3 ) + x 4 ] . 



(A4) 



Together with the moments, Table IB HB 3 1 these determine all tem- 
perature correction terms up to (9(0 3 ). 



distortion introduced by the scattering of photons off of electrons 
inside a moving cluster it is enough to have some way of computing 
the derivatives up to the desired order. Looking at the main scaling 
at high frequencies one can show that fl'/ipi ~ (— l)*e~ v /(l —e ~ x ) k+l . 
In general, one can also write 



(A5) 



where V k {x) is a polynomial of order k. 

With the Ansatz P k (x) = j} m =o aJ?*" 1 it is straightforward to 
find the recursion relation 



r !«=( m+ i) T r + (t- m )£» 



for the coefficient y,„ . The solution of this recursion is 



(A6) 



m I 



These coefficients are also known as Eulerian numbers, which 
can be generated in a very fast and numerically stable way using 
Eq. l |A6t . The derivatives of the Planckian spectrum are given by 



(-x)V v k A I 



(A8) 



for any k. Defining f H k = (- l)*x*e x /(l - e x ) k+l one can also write 
the derivatives as 



(A9) 



xd x n P \ 


= «i, 








2a2 
x a x n P \ 


= W 2 


[1 + e-*] 






3 o3 

x o\n P \ 


= <H 3 


1 +4<T A 


+ e~ 2x ] 




4 <j4 

x o x ti P \ 


= <H 4 


1 + lie" 


x + lle- 2 * 


-"I 


5 o5 

xro/i P \ 


= <H S 


1 +26e- 


x + 66<T 2r 


+ 26e~ 3x + e- 4x ] 


x 6 d 6 x n PI 


= K 6 


1 + 57e" 


* + 302e~ 2 


x + 302e~ 3 * + 51 e 



As these expressions show, for x » 1 the limiting behaviour is 
given by x*<9*«pi ~ "Hk ~ {—\)xre~ x . Similarly, in the limit x <k 1 
one has x k d k x n PI ~ (-l)*x*. For numerical applications these ex- 
pressions are very useful and easy to code up. 

To evaluate the sum Y m = ZtiT ' a l * JC*dJ«pi, which appears 
after collecting terms of different orders in e from the moments, 
I j , we can use the recursion relation 

/^ = T^[^^<0+/£i] ( A1 °) 

with /8i"j,', + 1)+1 = and K,„(x) = p[ m> (x) n n (x). This groups terms 
of similar order and yields numerically stable results, even if the 
coefficients at become very large (see Table IB it . 



A2 Closed form for the derivatives of a Planck spectrum 

In Eq. dA4t we give explicit expressions for the derivatives of the 
Planck spectrum up to the 8 th derivative. For fully analytic manipu- 
lations of the Boltzmann equation and when calculating the spectral 



APPENDIX B: AZIMUTHALLY AVERAGED CROSS 
SECTIONS 

The 5-dimensional collision integral can be reduced by 2 dimen- 
sions when considering the azimuthal symmetry of the scattering 
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process. In particular, we encounter the integrals 
d 2 cr 



/dcr A(p' di 
dO 7 An 



3o~T „? i 3 



1 - £ (a x ) + (al 



d 1 cr l 
d/i d/i' 



d _ o" 2 
d/jdyu' 



do- dtf>'d(f> 
dfl 7 An 



■ fifi 



, d cr 3o" T 3 , 3 r , , 



d;U d;U' 8 

f dcr d 



4/r 



(Bla) 



(Bib) 



(Blc) 



p 20"y" )-r-T7 + — rr - £V <r [P 2 (wO - 1] gin,y- ) 

' 16 



dyU dyU' 

9ffT „ 2 2 3 

16 s ' 



l+f|2+-fg-3<or 1B >[l-f| 



with g(ji,jj.') = ^(1 -/j 2 )(1 —fx' 1 ), £ = v'/(vy 2 K 2 ), and the averaged 
quantities (a sc ) = J(l - /j sc ) d(f>' dip/ An 2 = 1 - nn' and (a 2 c ) = 
/(l - A/ sc ) 2 d0/4yr 2 = <a- sc > 2 + 



APPENDIX C: MOMENT FOR MONOPOLE, DIPOLE 
AND QUADRUPOLE SCATTERING 

We are interested in kinematic corrections up to second order 
of the peculiar velocity. To this order, the CMB spectrum has a 
motion-induced dipole and quadrupole component inside the mov- 
ing frame. In n th order of 8 C we therefore need all moments /£., V\ . 
and K for s < 2n + 2. The required integrations were performed 
with the symbolic algebra package Mathematica, defining appro- 
priate replacement rules for all the appearing integrals in the Taylor 
expansion for small temperature. Note that for the Taylor expan- 
sion of the relativistic Maxwell-Boltzmann distribution function, 
Eq. ©, it is useful to transform to the variable £ = rj 2 /28 e . 

The moments related to the monopole part of the photon dis- 
tribution are summarized in Table IBTI for up to 10 th order temper- 
ature corrections, or in total 11 th order in temperature. Up to 5 th 
or der in temperatur e these expressions agree with the results given 
bv lltoh et alj dl998h when neglecting all recoil terms. 

The Tables IB2I and |B3| give the moments for the scattering 
of the dipolar and quadmpolar spectral anisotropy in the radiation 
field, but only up to for illustration. In SZpack we included 

tables up to 10 th order, similar to the monopole scattering case. 
The mo ments up to t he firs t order in temperature were also de- 
rived in IChluba et al.l d2012h . However, higher moments have not 
been given in the literature before. In this work, we use then to 
write down all kinematic correction terms up to second order of the 
peculiar velocity using explicit Lorentz transformations. 
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Table Bl. Moments for monopole scattering. Blank means the coefficient is zero. Within each row, the temperature order increases from left (oc 8 C ) to right 
(oc while in each column the order of the derivative of the blackbody spectrum increases from top (xd x np\) to bottom (x 22 3j 2 «pi). Each row therefore 
gives the coefficients a^"' for the moments i*j = 6 e 2JJ,_ ifi* 0£Yoo, while each column defines the spectral function Y„ = 2f"| 2 a^' ) x k d x iioo(x)Yoo- We did 
not list the moment /$„ = 7oo = ai "loo for Thomson scattering of the monopole in the table. 
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Table B2. Moments for dipole scattering. The table is organized similar to Table IB II Each row gives the coefficients d! m ' for the moments ti = 
#e 2" = o ^ °iYi m , while each column defines the spectral function D„ = 2|3) 2 JC*flJfli (x), with n\ (x) = £,„ n\ m (x)Y\ m . Since the dipole does not couple 
directly to the Thomson scattering cross section, all moments are at least first order in temperature. 
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Table B3. Moments for quadrapole scattering. The table is organized similar to Table IBT1 but we also add the coefficient ' = 1/10 to the table, which 
arises from Thomson scattering of the quadrupole anisotropy and is independent of temperature. Each row gives the moment coefficients for = 
2" = _i q^SiY^m, while each column defines the spectral function Q„ = Zj'Sf if' ^dfyvilx), with n 2 (x) = 2^ = _ 2 n 2m(x)Y 2m . 
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